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Abstract 

A random graph process, Gi{n), is a sequence of graphs on n vertices which begins with the 
edgeless graph, and where at each step a single edge is added according to a uniform distribution 
on the missing edges. It is well known that in such a process a giant component (of linear size) 
typically emerges after (l + o(l))^ edges (a phenomenon known as "the double jump"), i.e., at 
time t = 1 when using a timescale of n/2 edges in each step. 

We consider a generalization of this process, 5k (n), proposed by Itai Benjamini in order 
to model the spreading of an epidemic. This generalized process gives a weight of size 1 to 
missing edges between pairs of isolated vertices, and a weight of size K £ [0, oo) otherwise. This 
corresponds to a case where links are added between n initially isolated settlements, where the 
probability of a new link in each step is biased according to whether or not its two endpoint 
settlements are still isolated. 

Combining methods of [13j with analytical techniques, we describe the typical emerging time 
of a giant component in this process, tc{K), as the singularity point of a solution to a set of 
differential equations. We proceed to analyze these differential equations and obtain properties 
of Gk, and in particular, we show that tc{K) strictly decreases from | to as if increases 
from to oo, and that tc{K) — (1 + o(l)), where the o(l)-term tends to as iiT ^ oo. 
Numerical approximations of the differential equations agree both with computer simulations of 
the process QK{n) and with the analytical results. 

1 Introduction 

1.1 The Achlioptas problem and the biased process 

The random graph process on n vertices, Qi = Qi{n), introduced by Erdos and Renyi, is a sequence 
of (2) + 1 graphs, (^i'', . . . ,^1^2))^ where the Qi* is the edgeless graph on n vertices, and Qi"^ is 
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obtained by adding a random edge to ^i™ , chosen uniformly over all missing edges. 

A classical result of Erdos and Renyi ([9|) states that if T = and C < 1, then typically 
every connected component of Gi^ is of size O(logn), and if C > 1 then typically there is a single 
giant component of size 0(n) and every other component is of size O(logn). Thus, there is a phase 
transition (the "double jump") after (l+o(l))n/2 edges. For further discussion of this phenomenon, 
see, e.g., [2]. 

A well known problem, introduced by Achlioptas, discusses a scenario where two randomly 
chosen edges are presented at each step, out of which a single edge is chosen by some algorithm 
A. The goal of the algorithm is to postpone the emerging time of the giant component as much as 
possible. This was first examined by Bohman and Frieze in [3]. For additional results both on this 
problem (as an off-line and as an online problem) , as well as on the converse problem of creating a 
giant component ahead of time, see [1] , [5] , [6] , [7] . In [13], the authors describe a generic approach 
to analyzing the performance of algorithms for the mentioned Achlioptas problem. After applying 
Wormald's differential equation method for graph processes [16], the emerging time of the giant 
component is expressed as a singularity point to a differential equation. Using this method, the 
authors are able to provide bounds for the performance of several algorithms. 

In this paper, we study a natural generalization of the Erdos-Renyi random graph process, 
proposed by Itai Benjamini in order to model the spreading of an epidemic. This process of 
Benjamini lets "infected" clusters have either a larger or a smaller probability of increasing their 
size, depending on the value of an external continuous parameter K. This is achieved by embedding 
basic degree information (namely, whether or not a site is currently isolated) into the probability 
distribution over the missing edges. Note that the given model is the most natural of its kind with 
respect to the dependency of this probability distribution on the degrees of the vertices. Results on 
other, more complicated, models may be obtained via methods similar to the ones presented here. 

We combine methods of [13j together with analytical methods, with the same motivation of de- 
termining the critical time in which a giant component emerges in this process. As we mentioned, 
the generalized process we consider is a parameterized version, dependent on some K G [0, oo), 
which modifies the probability of each edge according to whether or not its endpoints are already 
connected. In the original Achlioptas problem, different algorithms can postpone the phase tran- 
sition or create it ahead of time, where the biased choice at each step is between precisely two 
randomly chosen edges. The process we study considers all missing edges when making its biased 
choice, and the phase transition is presented as a function of the continuous parameter K. 

By applying the powerful differential equation method, we were able to derive properties of 
our model directly from the system of coupled non- linear ordinary differential equations (ODEs). 
While some of these properties can be proved by relatively simple combinatorial arguments, cal- 
culating the precise asymptotic behavior of tdK) for K ^ 1 via combinatorial arguments seems 
challenging. Indeed, while at first glance this result appears as though it can be obtained by prob- 
abilistic arguments (e.g., using monotone coupling to variants of the model, tracking the structure 
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of the connected components throughout certain time intervals, etc.), we are not aware of any such 
derivation at the present time. In particular, the only way to obtain the 4/\/3-fC behavior of the 
blowup time is through the careful asymptotic analysis of the derived coupled differential equations 
(see, e.g., [1] for more on asymptotic analysis of ODEs). 

The generalized process can be efficiently implemented. In order to efficiently randomize the 
next edge, one needs to maintain the sets of isolated and non-isolated vertices, along with the set 
of edges already chosen. Our implementation runs in time 0(n log n) and requires 0{n) memory. 
The computer simulations show an excellent agreement with the numerical solutions of the ODEs 
on the one hand, and with the analytical results concerning tc{K) on the other hand. 

1.2 Notations and main results 

The biased graph process on n vertices, Qk, is the following generalization of the random graph 
process: as before, the initial graph Qk^ is the edgeless graph on n vertices, and Gk"^ is obtained 
by adding a single edge to Gk"^~^- The newly added edge is selected according to the following 
distribution on the missing edges: each edge between two isolated vertices is assigned a weight of 
1, and the weight K is assigned to all the remaining edges. Once there are less than 2 isolated 
vertices, the distribution on the missing edges is uniform. We extend the definition of Gk"^ to 
m > (2) by setting Gk"^ = -^n ^or every such m, where denotes the complete graph on n 
vertices. Furthermore, we use the notation Gk'^\h, where H is a graph on n vertices, to denote the 
biased graph process after T steps, starting from the initial graph H instead of the edgeless graph. 

Let Gxit) denote the biased process after scaling its time line by a factor of n/2: 

Since a choice of -R' = 1 is equivalent to the Erdos-Renyi random graph process Gi, the appearing 
time of the giant component in Gi is typically at t = 1. We study the effect that modifying K has 
on this critical point, tc{K), keeping in mind that, intuitively, decreasing the value of K should 
postpone the emerging point of a giant component and vice versa. 

Throughout the paper, we say that a random graph on n vertices satisfies some property with 
high probability, or almost surely, or that almost every graph process on n vertices satisfies a 
property, if the corresponding event has a probability which tends to 1 as n tends to infinity. 

Let G = {V,E) be a graph on |F| = n vertices. Let C = C{G) denote the set of connected 
components of G, and let Cy = Cy{G) denote the connected component of a given vertex v G V. 
We denote by Cq the set of isolated vertices in G, and by I{G) the fraction of isolated vertices: 

I{G) = ^ . 
n 

Notice that, if the graph after T steps (that is, once T edges have been added) contains |Co| isolated 
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vertices, then the probabihty that the next edge is a specific one between two isolated vertices is 



and the probabihty that it is a specific new edge between two vertices that are not both isolated is 
K times that value. 

The susceptibility of G, S{G), is defined to be the expected size of a connected component of a 
uniformly chosen vertex v £ V: 

5(G) = i 5] ICI = i 

v& C(LC(G) 

The susceptibility of G is closely coupled with the existence of a giant component in G. Indeed, the 
existence of a giant component of size an for some a > implies that S{G) > a^n, and conversely, 
if S{G) = i7(n) then at least one connected component is of linear size. Following the ideas of |13j , 
we characterize the behavior of the susceptibility along Qk, and in the process obtain the required 
results on tc{K). In order to do so, we need to examine the typical behavior of the number of 
isolated vertices along Qk- 

The fraction of isolated vertices in Qk-, I{QK{t)) has a value of 1 at t = 0, and decreases to 
over time. The following theorem summarizes the behavior of I{QK{t)): 

Theorem 1.1. For every K > Q and every C > 0, almost every biased process on n vertices Qk 
satisfies \I{QK{t)) -y{t)\ = 0{n~^l^) for every < t < G, where y is the solution to the differential 
equation: 

{i-y)K 

y^ + {l-y^)K ' . (1) 
y(0) = 1 

In the special case K = {), the above holds for C = 1. 

The following proposition analyzes the differential equation whose solution will prove to be a 
good approximation of the susceptibility along the biased process: 

Proposition 1.2. Let z = z{t) denote the solution for the following differential equation: 

= y2 + il-y2)K +1 , (2) 

z(0) = 1 

where y is the solution to the differential equation ([T|). For every > there exists a singularity 
point tc = tc{K) > 0, such that z{t) is continuous on [0,tc), and lim^^^- z{t) = oo. Furthermore, 
there exists some constant M > 0, independent of K , such that tc{K) G (0,M) for every K > Q. 

The next theorem implies that the singularity point of z{t), tc = tc{K), is the typical time at 
which a giant component emerges in a biased process: 



y' = y^-y>^^ 1 
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Theorem 1.3. For every e > and every K > 0, almost every biased process on n vertices Qk 
satisfies the following: 



1. Subcritical phase: \S{QKit)) — z{t)\ = o(l) for every t G [0, — e]; where z{t) is the solution 
for the differential equation ([2]) , and tc is its singularity point as defined in Proposition \1.SX 
Furthermore, for every t E [0,tc — £]> ihe largest component ofGxit) is of size O(logn). 

2. Supercritical phase: S{QK{tc + ^)) = ^{n). 

Altogether, the appearance of a giant component in Qk is almost surely at time tc- In the special 

{ l + t ift<l 

case K = f), the above holds when replacing z[t) with the function z{t) = < . 

I 372^1 ^/*>1 

Finally, the behavior of tc{K), the typical point of the phase transition in is characterized 
by the following theorem: 

Theorem 1.4. Let tc{K) denote the singularity point of the solution to the ODE as defined 
in Proposition \l.^ Then tc{K) is continuous and strictly monotone decreasing as a function of K 
and satisfies: 

tM = I 

where the o{\)-term tends to as K ^ oo. 

The rest of the paper is organized as follows: in Section [2] we study general properties of the dis- 
tribution of Qk- Namely, we study the relation between Qk and Qi, and describe an approximated 
process, Qk, which is easier to analyze. 

In Section [3] we prove Theorem 1 1 . 1 1 and analyze the solution to equation ([1]), which characterizes 
the behavior of the isolation ratio throughout the biased process. Prop osit ion 11.21 and Theorem 11.31 
are both proved in Section where computer simulations are also included. 

In Section [5] we analyze the asymptotic behavior of tc{K) for large values of K, and prove 
Theorem 11.41 The final section [U] is devoted to open problems and concluding remarks. 



2 Dominating and approximate distributions for 

2.1 The relation between Qk and Qi 

A main tool in understanding and analyzing the biased process Qk is the stochastic domination re- 
lation between Qk and time-stretched versions of the Erdos-Renyi process Qi- In order to formalize 
and prove this relation, we consider a wider family of graph processes, defined as follows: 
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Definition. Let M G N. An M -bounded weighted graph process on n vertices, 7i = 'H{n), is 
an infinite sequence of graphs on n vertices, {7i^,Ti.^, . . .), where Ti^ is some fixed initial graph, and 
TC" is generated from hy adding one edge at random, as follows: the probability of adding the 

edge e to H^~^ is proportional to some weight function Wt{e), satisfying: 

max Wt(e) < M min Wt(e) . 

If for some > Ti^ = Kn, we define = TC^ = Kn for every t > v. 

Clearly, the biased process Qk is an M-bounded weigiited grapli process which starts from the 
edgeless graph, where M = [max{i(', -^j] . We are interested in the relation between the probability 
that Qi satisfies some monotone graph property A (a set of graphs closed under isomorphism and 
under the addition of edges) and the corresponding probability of Qk- The following theorem 
formalizes the stochastic domination of the original graph process Qi on M-bounded weighted 
graph processes: 

Theorem 2.1. Let Ti. denote an M-bounded weighted graph process on n vertices, and let A denote 
a monotone increasing property of graphs on n vertices. The following statements hold for any 
t e N; 

Pr[7^* e ^] < V^[gi^''%. G A] , (3) 
Pr[gi%o G ^] < Pr[?^^^* G A] . (4) 

We need the following lemma, which was first proved in [T3] in a slightly different setting. For 
the sake of completeness, we include a short proof of the lemma using the Max-Flow-Min-Cut 
Theorem (a relation which was first observed in |12j): 

Lemma 2.2. Let U, V be two finite sets, and let R C U x V denote a relation on U, V . Let /i and 
V denote probability measures on U and V respectively, such that the following inequality holds for 
every A (Z U: 

fJ-{A) < I'dy G V : xRy for some x G A}) . (5) 

Then there exists a coupling ip of fijU whose support is contained in R. That is, there is a joint 
distribution if on U,V satisfying the following two properties: 

1. The marginal distributions of on A and B are fi and v respectively. 

2. For every {x,y) G U x V such that -^{xRy), ip{x,y) = 0. 

Furthermore, the distribution ip can be found in time polynomial in\U\ + \V\. 

Proof. Consider a directed weighted graph G on the vertex set C/ U y U {s, F[, with the following 
set of edges and capacities (see Figure [1]) : 



6 



x'Ry' 



Figure 1: The network graph G in Lemma |2.2[ 

1. For every x £ U and y £ V such that xRy, place the edge {x, y) with a capacity of oo. 

2. For every x £ U, place the edge {s,x) with a capacity of 

3. For every y £V, place the edge {y,t) with a capacity of 

Note that a flow of 1 from s to t in the network defined above provides the desired distribution ip, 
by defining ip{x,y) = f{x,y) where / is such a flow. The flow value of 1 implies the first property 
required from ip. The second property follows from the fact that, if ^{xRy) for some x £ U and 
y £ V, then the edge {x,y) is not in G. Altogether, if indeed G has a flow of 1, then a joint 
distribution ip satisfying the above properties can be computed in time polynomial in \U\ + \ V\. 

The fact G has a flow of 1 follows essentially from the proof of Hall's Theorem using the Max- 
Flow-Min-Cut Theorem (see, e.g., |10]). where inequality ([5]) replaces Hall's criteria for a maximal 
matching. Let f{G) denote the maximal flow from s to t in G. By the Max-Flow- Min-Cut Theorem, 
/ corresponds to a minimal (S, T) cut, where s £ S and t £ T. Hence, if we define A = S CiU and 
B = S nV, the choice of the capacity oo for the edges between U, V implies that N{A) nV C B, 
where A^(^) is the set of neighbors of A in the graph G. Hence: 

f = fi{U\A) + v{B) > ii{U \ A) + v{N{A) r\V)> fiiU \ A) + fi{A) = 1 , 

where the last inequality is by the assumption ([5]). ■ 

Proof of Theorem \2.1\ Let Q denote two graph processes starting from the same initial graph. 
We wish to prove that, under certain conditions, the following inequality holds for every t S N: 

Pr[P* G ^] < Pr[Q^^* £ A] . (6) 

The proof will follow from a coupling of the two processes V and Q, such that for every instance of 
the coupling and every t G N, is a subset of Q**^. More precisely, we define a joint distribution 
on the processes ("P, Q) in the following manner: at each step t G N, we generate according to 
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its usual law, then generate the pairs 

g(t_i)M+i^ . . . , Q**^ according to some distribution, such that 
E{V*) C E{Q^^^) and the marginal distribution of Q{t-^)M+i ^ _ _ _ ^ QtM jg correct one. Clearly, 
such a construction immediately proves ([6]). 

For every t gN and i £ {!,..., M}, let XI and It denote the sets of missing edges of q(*-i)a^+*^i 
and respectively. Let ej G denote the t-th edge added to V, and let G X^* denote the 
((t - 1)M + i)-th edge added to Q. 

Recalling that, by definition, "P*^ = Q*^, let t > 1 and assume that we have constructed the 
above coupling up to the point Q(*~^)^^). In particular, E{V*~^) C E{Q^^~^^^), and hence 

Xl C Yt. 

Note that in the special case M > \X^\, in which Q**''^ is by definition the complete graph, we 
can clearly generate q(*~^)*^+^, . . . , q(^-'^)^+\^\ according to the usual law of Q and maintain the 
properties required from the coupling. Assume therefore that M < 

We now wish to extend the coupling to ("P*^'^, Q*) using the joint distribution (p of Lemma |2.2[ 
To do so, take: 

' U = Yt 

V = {XDm 

/u(e) = Pr[ej = e] for every e ^Yt , (7) 

v{F) = Pr[(/i,...,/AO for every F G (X1)m 

R = {(e, F)eYtX {XDm ■ e e {Yt\ X}) or e G F} 

where {XI)m denotes the set of ordered subsets of M edges of X} . 

If inequality ([5]) holds for every A <ZU , then (p associates every edge e G It with a distribution 
on Q(*- 1)^^+1,..., Q**^ which ensures that e G Q**^, by the second property of the lemma. The 
choice of v, along with the first property of the lemma, ensures that Q{*-i)A^+i, . . . , QiM ^iii have 
the correct marginal distributions. 

Clearly, ([5]) holds for every A C U such that Ari{Yt\ X^) ^ 0, since by the definition of R, the 
right hand side of ([5]) equals 1 in this case. It remains to show that inequality ([5]) indeed holds for 
every A C Xf when V and Q play the roles of the processes ^i|-ho and 7i. This is proved in the 
following two claims: 

Claim 2.3. Let V = Ti. and Q = Gilno, and define U,V, fijU, R as in ([7]). Then inequality ([5]) 
holds for every A C U. 

Proof. To simplify the notations, let X = X} and let Y = Yt. By the above choices for V and Q 
and by the definition of the M-bounded weighted process, we have: 

v{F) = Pr[(/i, . . . , /f ) = F] = — ^ for every F G {X)m , 

WAe) 

//(e) = Pr[et = e] = = „. , for every e G F , (8) 
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where, without loss of generahty: 

minW((e) = 1 , maxWt{e) < M . 

Substituting the values of ^, v, inequality ^ takes the form: for every A C X, 

(\X\-\A\\ 

Pr[e, eA]<l- ^-M^ . (9) 

\Al) 

Take AG X, and assume that A^9,X (otherwise Q trivially holds). By ([8]), we have: 

Pr[e. E A] - ^Fr[e, - e] - ^-^^ < |y| + (m - 1^1 ' 

where the last inequality is by the fact that 1 < Wt{e) < M for every e G y, along with the 
inequality > ^ for every < a < 6 and c > 0. Recalling that X C y, we get: 

^'1'' ^1 ^ |X|+'(M-1)|A| ■ <1°> 



(11) 



Combining pOj) with the fact that: 
we obtain the following sufficient condition for ([9]): 

- |x| - |A| + M|^| ■ ^ ' 



Set /3 = 1 — and recall that < /3 < 1. Inequality (fT2|) takes the following form: 



dM ^ P 



p+{l-(5)M ' 
or equivalently: 

P^-M _ (1 _ M)/3 - Af > . (13) 

Defining g{x) = x^~^^ - (1 - M)x - M, it is easy to verify that g'{x) < for < x < 1 provided 
that M > 1, and that g{l) = 0. Hence, inequality (fTSll indeed holds, and inequality ([9]) follows, as 
required. ■ 

Claim 2.4. Let V = ^iI-ho and Q = 7i, and define U,V, fijU, R as in ([7]). Then inequality ([5]) 
holds for every A C U. 

Proof. Following the notation of the previous claim, define X = Xl and Y = Yt. As /x is the 
uniform distribution on Y , inequality ([5]) takes the following form: for every A C X, 

|^<l-Pr[//,...,/,^^^] . (14) 
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Let A C X, and we may again assume A ^ (}>,X otherwise (|14p trivially holds. For i = 1, . . . , M, 
let denote the weight function on the set of edges XI. According to this notation, 

the right hand side of (jl4p satisfies: 



M 



l-Fr[f^,...jr^A] = l-llFT[fl^A\f^,...Ji-'^A] = 

i=l 

By the same argument used in Claim we reduce the expression ^^^^ — * ^ \ by assigning the 
value 1 to the weights of A and M to the rest. Hence, 

i ^A\>i iii^i \xi\M + \A\{i-M)) Vi^i^-i^i^ + i^iy ■ 

Recalling that X C Y and setting /? = 1 — (0 < /3 < 1)) the claim then follows from the next 
inequality: 

which holds for every /3 > whenever M > 1. ■ 

Combining Claim [231 with the arguments preceding it and Lemma 12.21 yields that it is possible 
to extend the coupling to (7^*^, Q*) when V = Tl and Q = GilnOj thereby completing the induction 
argument. Therefore, ^ holds for every t G N. Similarly, combining Claim 12.41 with the above 
arguments yields that dH) holds for every t G N. This completes the proof of the theorem. ■ 

Remark 2.5: Theorem [2T] required that one of the processes, 7i, chooses each edge according to 
an A/-bounded distribution, whereas the other chooses each edge according to a uniform distribu- 
tion. It is not difficult to construct an example showing that this requirement cannot be replaced 
by the condition, that the maximal ratio between the weights of the two processes at each step is 
at most M. 

Applying Theorem 12. II on the biased graph process Gk gives the following immediate corollary: 

Corollary 2.6. Let A denote a monotone increasing property of graphs on n vertices, and let t > 0. 
Then for every K > 0, the following two statements hold: 

1- IfGK{t) almost surely satisfies A, then Qi{\max{K, ^}~\t) almost surely satisfies A. 
2. If Qi{t) almost surely satisfies A, then [maxjii', -^Hi) almost surely satisfies A. 
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2.2 The approximate biased graph process 



In order to simplify the proofs of Theorems 11.11 and 11.31 we consider a variant of the biased process 
on n vertices, Qk, which we dub an approximate biased process. At each step, a random ordered 
pair of vertices {u, v) £ (where V is the set of vertices) is chosen out of the possible pairs, 
according to the following distribution: a pair of isolated vertices has a weight of 1, whereas all 
other pairs have weights of K. If the chosen pair corresponds to a (self) loop or to an edge which 
already exists in Qxi no edge is added in this step. 

The following claim implies that it is sufficient to prove Theorems 11.11 and 11.31 for the approximate 
model: 

Claim 2.7. Let C > 0, and let {At '■ < t < C} denote a family of properties of graphs on n 
vertices. If for almost every approximate biased process Qk, Gxit) satisfies At for every < t <T, 
then for almost every biased process Qk, Gxit) satisfies At for every < t <T as well. 

Proof. Let Bj denote the event that the ordered pair, chosen in the j-th step, was not added to 
Gk, either being a loop or already belonging to Gk- Fix C > 0, and let T = C^. The probability 
of the event Bj (1 < j < T) satisfies: 

^ r„ 1 1 2(j/- - 1) C+l 

and, as these events are independent, we apply the well known bound 1 — x > e~^/(^~^^ for 
< X < 1, and obtain: 

Pr[AjliB,] >[l- -^j > exp ( -^3^ I ^ exp(-C(C+ 1)) . 

. — I 

Notice that if we condition on the event Aj^^Bi, then the two graph sequences {Gk ,---,Gk ) 

and {Gk^, ■ ■ ■ ,Gk^) have the same joint distribution. Since there is a fixed lower bound on the 

— 1 — T 

probability for this event, any statement on [Gk , ■ ■ ■ , Gk ) which holds almost surely, also holds 
almost surely for {Gk^, • • • , Gk^)- The result follows. ■ 

Remark 2.8: In the above claim we used the low probability for a step to get omitted in order 
to show that statements that hold almost surely for Gk can be derived from such results on Gk- 
However, it is worth noting that the two processes are much closer than that; if we condition that 
Gk has M edges, then clearly it is distributed as Gk - It is not difficult to show that, by the low 
probability for omitting a step, GKit) has the same distribution as Gk ((1 + o{l))t). 
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3 The behavior of the isolation ratio 



3.1 Proof of Theorem fTTl 

We begin with the special case K = 0, which we prove directly on Qk- In this case, the differential 
equation ([T|) takes the simple form: 

y' = -l, y(0) = 1 , 

and hence its unique solution is y{t) = 1 — t. Notice that, as K equals 0, the biased process connects 
two isolated vertices at each step with probability 1, as long as two such vertices exist. Hence, at 
time < t < 1, there are [tn/2\ edges which are vertex disjoint in pairs. Thus the ratio of isolated 
vertices, /(^^^(t)) = \CQ{QK{t))\/n, satisfies: 

\i{g^{t))-{i-t)\<^ , 

n 

where the ^-term is the rounding error. 

We are left with the case -fC > 0. Fix C > 0, and let G ~ Qxit) denote a graph at some point 
t < C along the approximate biased process Qk- We examine the effect that a single step of Gk 
has on the ratio isolated vertices, I{G) = \CQ{G)\/n. 

Set / = I{G); the total of the weights assigned to all ordered pairs is (/^ + K{1 — 
Hence, with probability j2_^j^(^i"_f2-^ the chosen edge is between two formerly isolated vertices, and 
with probability precisely one end point of the chosen edge was formerly isolated. Let 

G' denote the graph after performing the above step. The expected change in the isolation ratio 
between G and G' satisfies: 



n 



2 / (1 - I)K I/n , 

'ny-' P + K{l-P)~ P + K{l-P)' ■ ^^^^ 



Define: 



erry(/,K) = j^^^^^f^^, (17) 

and notice that the denominator of erry(/, K) lies between 1 and K for every value of I. Hence, as 
K>0: 

I In 

As G' corresponds to QK{t + ;^)) we choose At = 2/n and rewrite ()16p in the following form: 

/(fe(t + At))-/(fe(t)) ^ (1 - I)K 

At p + K{l-P) -^ + «"y^^'^^- ^^^^ 

Notice that the left hand side of (fT9ll resembles suggesting that the expected change in the 
isolation ratio I{QK{t)) is linked with the solution to the differential equation ([1]). Notice that y' 
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is C°°, and hence there is a unique solution to ([T|): we analyze this solution in the next subsection, 
and now turn to show that indeed it approximates I{Qi^{t)). This will follow from a general 
purpose theorem of [16] (Theorem 5.1), which we reformulate according to our needs, for the sake 
of simplicity (a simpler version of the theorem appears in [l5]). If y = {yi, . . . ,yk), we use the 
notation (x, cy) to describe the tuple (x, cyi, . . . , cy^)- 

Theorem 3.1 ([H]). Let Y = Y{n) = (Yi, . . . , 1/) denote I > 1 functions from graphs on n vertices 
to the real interval [—Con, Con], where Co > is some constant, and let f = (/i, ■ ■ ■ , fi) denote I 
functions in m!''^'^ M. Let {H''"} denote a graph process on n vertices beginning with the edgeless 
graph. Let T) C M'^"*^ he a bounded connected open set such that (0, — ^^-^) £ T^, and let Tx>{H) 
denote the minimal time T such that (— , ' ) no longer belongs to V. Assume the following: 

1. (Boundedness Hypothesis) For some function /? = /3(n) > 1, \Y{H'^^'^) — Y{H'^)\ao < j3 for 
every T < Tjj^H). 

2. (Trend Hypothesis) For some function Ai = Ai(n) = o(l), the following holds for every 



T<Td{H): E(Y{H^+')-YiH^))-fi^,^ 



< Ai 



1 

oo 



3. (Lipschitz Hypothesis) There exists a constant L > such that the following holds for every 
x,ye V: \ f{x) - /(y)|oo < L|x - yloo- 



Then the following holds: 



1. There exists a unique solution to the following system of I differential equations: 

^ = /i(x,n), ^ = /Kx,n) 

where u denotes (ui, . . . ,ui). Let ui{x), . . . ,ui{x) denote this solution. 

2. Let A > Ai, A = o(l), and let a > be such that {x,u(x)) is at least CiA away from the 
boundary of T> for every x £ [0,cr], where Ci > is some constant. Then there exists a 
constant C2 > 0, such that with probability 1 — 0{j exp(— ^^)), the following holds for every 

0<T <an and every l<i<l: |^^^^ - ^ C'aA. 

Remark 3.2: Note that for the purpose of proving Theorem 11.11 we will only need the special 
case ^ = 1 of Theorem 13.11 above. However, in Sectional when analyzing the critical point for the 
emerging time of the giant component, we apply the above theorem for / = 2. 

We claim that the following substitution completes the proof of Theorem 11.11 

Y{G) = \Co{G)\ , 

r + (1 - y^)K 
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For the set V we choose a bounded connected open set containing the rectangle [0, C] x [0,1]. 

— T 

Notice that < < 1 for every graph G on n vertices, hence (^, ^ ) belongs to V for every 

1 < T < Cn/2 and Td{Qk) is at least C. 

— T+l — T 

Indeed, \Y{Qk ) — Y{Qk )| < 2, thus a choice of /3 = 2 confirms the Boundedness Hypothesis. 
Next, (HI]), (HZ]) and ([HD imply that 

IE (F(eV) - Y(g'/)) - /(^, = 2e„,(B^, = 0(l/„) . 

V J n n n 

Thus, setting Ai = ;^ verifies the Trend Hypothesis. Finally, f{x,y) is clearly (recall that 
K > 0) and hence satisfies the Lipschitz condition. 

Proposition 13.3^ proved in the next subsection by analyzing the differential equation ([T]), states 
that the solution u{x) is bounded between exp(—x/Q) and exp(—Qx). Hence, for every < x < C 
the solution u{x) is at least a constant away from the boundary of and we can easily choose 
a = C (regardless of our choice of A). 

Altogether, a choice of A = n^^/^ implies that with probability 1 — 0(n^/^ exp(— n^/^)) = 
1 — o(exp(— in^/^), the quantity Y{G) satisfies: 

^^=^(^)+0(n-V4), 
n n 

for every < T < Cn/2, where u is the (unique) solution to the equation ^ = f{x,u). 

Scaling the time to units of n/2 edges, we obtain that, with probability 1 — o(exp(— ^n^/^), the 
following holds for every < t < C: 

\I{gK{t))-y{t)\=0{n~'/') , 

where y{x) = u{x/2). Hence, y is the unique solution to the equation: 

dy 1 du (1 — y)K ^ 

dx ~ 2d^ ~ y2 + (1 - y'^)K ~ ' 

with the starting condition y(0) = n(0) = 1, completing the proof. ■ 
3.2 The behavior of y(t) when K > 

The behavior of y{t) along the biased process is crucial to the understanding of how the susceptibility 
grows, as we show in the next section. The behavior of the ratio of isolated vertices (and the 
corresponding function y{t)) was already stated for the case K = 0. We thus assume K > 0, and 
set Q = \max{K, ^}~\- By Corollarv 12.61 we obtain that for every e > and every C > 0, the 
following holds with high probability for every t £ [0, C]: 

e-Q* -£< liGKit)) < e-*/^ + e . (20) 
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This follows from the uniform continuity of the functions exp(— t/Q) and exp{—Qt), which describe 
I{Qi{t/Q)) and I{Qi{Qt)) respectively. 

We next show that a stronger result than (j20p can be easily derived directly from the differential 
analysis of y{t), as stated by the next proposition: 

Proposition 3.3. The ODE ([T]) for > has the following properties: 

1. Its solutions are strictly monotone increasing in K. That is, if yi{t) , y2{t) are the solutions 
that correspond to Ki < K2, then yi{t) < y2{t) for all t > 0. 

2. The solution y{t) is strictly monotone decreasing in t, and satisfies: 

{ ri r^i l ^ ^ { JTITX] ' (21) 

\ mm|l,Aj/ \ maxjljAjy 

for all t > 0. Furthermore, the inequalities are strict whenever K ^ 1. 

Proof. We begin by proving an inequality analogous to ([20|) . stating that whenever K ^ 1, y{t) is 
strictly monotone decreasing in t and satisfies: 

exp(-Qt) < y{t) < exp{-t/Q) (22) 

for any t > 0, where Q = max{-^,i^}. Inequality ()2ip will follow directly from (j22l) once we prove 
that y is strictly monotone increasing in K, since exp(— t) is the solution to the ODE ([1]) for K = 1. 

Suppose K > 1 and let u{t) = exp{— t/Q}. First, we prove that u satisfies the inequality 

(l-u)K , , , 

li^ + (1 — li^jA 

1 U 

Indeed, u' = -— exp{-t/Q} = - — . Therefore, ([23]) holds iff 

u (1 — u)K 

-j^> sin 2^^ -!' forO<n<l. 24 
Q + (1 — n^JA 

which, after some manipulations, is equivalent to the inequality 

Pq{u) = {K - l)u^ - Q{K - l)u + K{Q - 1) > 0, for < u < 1. 

The quadratic polynomial Pq{u) satisfies Pq{u) +00 as u — > ±00, because K > 1. Furthermore, 
if Q > 1 then pq{0) = K{Q — 1) > and Pq{1) = Q — 1 > 0. Therefore, if Pq{u) does not 
have a root in the interval (0, 1) then the inequality holds. In the special case of Q = X we have 

Pk{u) = {K-l) (n^ - Ku + K). The roots oipxiu) are u± = ^ ±\^\J K"^ - iK. For 1 < K < 4 
there are no real roots. For K = 4 there is a double root at u = 2 and for K > 4 there are two 
distinct roots u± > 1. Therefore, the inequality (p3]) holds for all values of Q = > 1 as asserted. 
By the equivalent inequality (I24p (recall that < y < 1 for t > 0), we obtain that y'{t) < for 
every t >0, and hence y{t) is strictly monotone decreasing. 
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To complete the proof of the upper bound of (I22]) . let f(w) = — k 1. The 

u;^ + (1 — w'^jK 

differential equation ([1]) and inequality (j23p imply that the functions y{t) and u(t) satisfy for t > 

y' = f{y), 

u > f{u), 

together with the mutual initial condition y(0) = u{0) = 1. By standard analytical considerations, 
this implies that u{t) > y{t) for all i > (to see this, note that u{t) > y{t) for t £ [0, e] and 
some small e > 0, and that the smallest point t* > e satisfying u{t*) = y{t*) cannot satisfy 
u'{t*) - y'{t*) > f{u{t*)) - f{y{t*)) = 0, yielding a contradiction). 

The lower bound in (j22p and the case < K < 1 are similar. 

To prove that y is strictly monotone increasing in K, take < Ki < K2, and let yi and 2/2 the 
solutions of equation ([T|) for Ki and K2 respectively. Considering /, defined as above, as a function 
of both w and K, f satisfies: 

5/ ^ {l-w)w^ 

dK^ ' ' {w-^ + {I - w'^)KY ' 

hence ^{w,K) > for every < ?i; < 1. Recall that inequality (f22]) guarantees that y{t) > for 
every t. Therefore, for any t > 0, f(y{t),K) is monotone increasing in K, and in particular, for 
any t > we have y[ = g{yi) and y2 > 5(2/2), where g{y) = f{y,K). Thus, the above argument for 
proving inequality ([22|) completes the proof. ■ 

We note that the bounds of (j21l) are much weaker than the estimations which can be obtained 
by examining the asymptotic behavior of the ODE ([1]), as we proceed to do in Section [5l 



4 The susceptibility of the biased process 
4.1 Proof of Proposition [T72] 

We prove Proposition 11.21 by showing that there exists an M > such that < tc{K) < M for 
every K > 0. The fact that z{t) is continuous on [0,tc) follows from standard considerations in 
differential analysis. The case K = 1 is trivial, as in this case equation ([2]) takes the simple form 

z' = , z{0) = 1 , 

hence its solution is 1/(1 — t) and tc{l) = 1. 

Take K > 1, and let z{t) denote the (unique) solution to equation ([2]). Recall that by Proposition 
13.31 y{t), the solution to equation ([T]), satisfies: 

e-* < y{t) < e-*/^ 
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K 

for every t > 0. Thus, y{t) > for every t > 0, and we obtain that ^5 ij— — > 1 for every 

r + (1 - y^)K 

t > 0. Therefore: 

z' > , z(0) = 1 , 

and by the method used in the proof of Proposition 13. 3|. we obtain that z{t) > 1/(1 — t) for every 
t > 0. We deduce that ^ <tc{K)<l for every K > I. 

Let < X < 1, and again let y{t) and z{t) denote the solutions to equations ([T]) and ([2|) 
respectively. Recalling that 2/(0) = 1 and y{t) decreases to as t ^ 00, let t* > be such that 
y{t*) = y/K. It is easy to verify that the derivative of /(y) = ^a^^^^-y^jA" ~ ^ strictly negative for 
every < y < 1. Thus, ([1]) implies that for every < t < t*, 

y'(t) = /(y) < = ^-^^ ■ (25) 

By the Mean Value Theorem, we obtain that: 

f < = , ' , < 2 , (26) 

~ 1 \-^fK 1 3_ - ' ^ ' 

where the last inequality is by the fact that the function ^ — monotone increasing from 

^ to 00 as goes from to 1. In addition, by (l2|), z'(t) > for every t > (this applies to every 
IT > 0). Therefore: 

z{f) > z{0) = 1 . (27) 

Hence: 

. ^ K 1 

r + (i-r)-?^ 2 

for every t > t* . By ([27]) and defining z{t) = z{t + t*) gives: 

z' > , z{0) > 1 , (29) 

and by comparing z with the function we deduce that z has a singularity point tc ^ 2. Alto- 
gether, z{t) satisfies tc < tc + t* < 4, completing the proof. ■ 

Remark 4.1: Theorem 11.41 states that tc{K) is in fact monotone decreasing as a function of K 
and thus bounded by tc(0) = |. We note that Theorem 11.31 used in the proof of Theorem 11.4^ 
requires a bound on tc{K). 



4.2 The susceptibility at the subcritical phase 

In this subsection we combine Theorem 12.11 with methods from [13] to prove part [1] of Theorem ll.31 
The proof relies on the fact that a bounded susceptibility ensures a logarithmic upper bound for 
the components, and this in turn ensures that the solution to the differential equation ^ stays a 
good approximation for the susceptibility. Hence, the fact that S{GK{t)) is approximated by z{t) 
and the fact that every component is of size O(lnn) ensure each other along the subcritical phase. 
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As we next show, the special case K = again follows from the behavior of S{Qi). Indeed, ([2|) 
takes the following simple form when K = 0: 



z' = l, z{0) = 1 , 

and hence its unique solution is z(t) = 1 + z. By definition, as = 0, as long as there exists a pair 
of isolated vertices in the graph, such a pair forms the next edge. This amounts to trading two 
components of size 1 with one component of size 2, and hence at each such step the susceptibility, 
being the sum of squares of the component sizes, increases by ^. Thus, at time < t < 1, the 
susceptibility equals 1 + ^[tn/2\, hence: 

\sigK{t))-ii + t)\< - , 

n 

where the —term is the rounding error. As of this point, the process is equivalent to Gi on [^^^J 

vertices, [n/2j of which represent components of size 2. Hence, the result on S{G) from this point 

on is derived from the case = 1 of Theorem 11.31 The solution for the differential equation ([2]) 

when K equals 1 is z = hence, we obtain that the susceptibility is within o(l) of the function 
? — 2 1 r J. > 1 

Throughout the remainder of the proof, assume therefore that K > 0. We need the following 
definition: a graph G is said to have a k,c component tail if, for every s, the probability that a 
uniformly chosen vertex has a component of size at least s is at most /cexp(— cs): 

Pr [ICJ > s] < ke-"' . 
veV(G) 

Notice that if a graph G on n vertices has this property then its largest component is of size at 
most max{l, l/c}lnn, provided that n is sufficiently large. 

The following theorem ([13], Theorem 3.1) is crucial to the proof: 

Theorem 4.2 ([13J). Let £,k,c be positive real numbers. Let G be a graph on n vertices with a 
k, c component tail. Let H be a random graph with edge probability p = ^ on the same vertex set, 
where t is fixed. Set G^ = G U H. 

1. (Subcritical) Assume S{G) < i for all n. Let t£ < 1. Then there exist k~^ , c"*" (dependent 
on k,c,l,t but not on n nor G) such that G^ almost surely has a fc^,c^ component tail. In 
particular, all components have size O(lnn). 

2. (Supercritical) Assume S{G) > i. Let ti > 1. Then G^ almost surely has a giant component. 
More precisely, there exists 7 > (dependent on k,c,£,t but not on n nor G) such that G^ 
has a component of size at least ^n. 

The following lemma proves part [1] of Theorem 11.31 (analogous to Theorem 4.2 of jl3j): 

Lemma 4.3. Let e > 0. For every t £ [0,tc — e], where tc is as defined in Proposition the 
following two conditions hold almost surely: 
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1. (Susceptibility Approximation) \S{GKit)) — z{'t)\ = o(l), where z{t) is the solution to the 
differential equation 

2. (Small Components) Gxit) has a k,c component tail, for some k = k(t),c = c{t). 
Proof. Set Q = [max{-^, K}~\ , and let: 

L = \izitc -e)+ e)Qtc\ , t, = (j = 0, . . . , L) . (30) 

Assume by induction that properties [1] and [2] of the lemma almost surely hold for t £ [0, tj] 
for some < j < L (the case j = is trivial). We show that with probability 1 — o(l) properties 
[T] and [2] hold for the interval as well, and the lemma follows from a union bound on the 

complement events. 

First, consider the Small Components property. By the induction hypothesis and the mono- 
tonicity of z{t) (recall that z'{t) > for all t > 0, as mentioned in the proof of Proposition ll.2p . 

SiGKitj)) < Z{tj) + o(l) < z{tL) + 0(1) . 

Let: 

i = z{tL)+£ , i = Qj- . (31) 

With high probability, GK{tj) has a k,c component tail for some constants k,c (dependent on tj), 
SiGxitj)) < i and t£ < < 1. Thus, Theorem 14.21 implies that ^i(Olg^(tj) almost surely has 
a k~^,c~^ component tail for some constants k~^,c^. By Theorem [2?T] (inequality ([3])) we obtain 
that GK{t/Q)\gjf(tj) almost surely has a k~^,c'^ component tail, hence Gxitj-^-i) satisfies the Small 
Components property. Thus, in particular, with high probability all the components of GK{tj+i) 
are of size O(lnn). 

We next analyze the change in the susceptibility along the approximate biased process. Since 
Gk^ is equivalent to Gk for some T' > T (where T' — T is the number of redundant steps in 
Gk), the Small Components property which we proved above implies that, with high probability, 
the largest component of Gxitj+i) is of size at most a Inn for some a > 0. Assume therefore that 
this indeed is the case. 

Let G denote an instance of Gxit) on n vertices for some t < tj+i, and set / = I{G) and 
S = S{G). Recall that the total of the weights assigned to all ordered pairs is (/^ + K{1 — /^))n^, 
and that in case that the chosen pair in the next step is a loop or an edge which already exists 
in G, the step is omitted. Let G' denote the graph after the next step, and let Aci,C2 denote the 
event in which the newly chosen pair joins the two components Ci,C2 € C{G). If G' equals G or 
if the new edge is an internal edge of some component C G C(G) (the event Ac\c occurred), then 
S{G') = S{G). Otherwise, a new edge between Ci,C2 € C{G) implies: 

S{G') = - V \Cf = S + -\Ci\\C2\. (32) 
n ^-^ n 

ceCiG') 
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Hence: 

E(.(g) 5(G)) |C.r.|P.[Ac.,c.|. (33) 

' CidC C2€C 

There are two cases to consider, according to which we divide the summands of the right hand side 
of (133]) to two terms, Ai and A2. In the first case, a new edge (u, f) is added between two distinct 
and formerly isolated vertices u^v. Since |C„| = \Cy\ = 1, the contribution of this case to the sum 
in (|33p is its probability: 

Ai = ^ Pr[Ac„,cJ = j2 , Q _ j2^K ~ P'^fU-eCo^c^.cJ • (34) 

In remaining case, ^Ci,C2 occurs where Ci 7^ C2 and at least one of the components Ci,C2 is not 
an isolated vertex. The contribution to the right hand side of (j33p in this case is: 



' iH (/2 + (1 - P)K)v? /2 + (1 - P)K ^ (/2 + (1 - P)K)n? ■ ^ ' 

C2eC |C|>i 



Define: 



err,(I,E:) = Pr[U„gCo^C„.cJ + ^ (J2 + (f {2)^)^2 = 

|q>i 

= ^ yl^^ ^ (36) 

P + {l-P)K ^ n"^ P + {l-P)Kn' ^ ' 

Combining (f34l) . (f35]l and (f36]l . the following holds: 

+ = i2 + (i-P)K + /2 + (1 - P)K - ^) ' 

C2ec 

and therefore: 

E(5(G') - 5(G)) (1 - K)P K 



2/n P + {l-P)K P + {l-P)K 



S\G)-eTT,{I,K) . (37) 



This explains the choice of the differential equation ([2]), provided that errz is redundant; this is 
ensured by the fact that the largest component of G is of size at most a Inn. To see this, notice 
that G satisfies ^ YliC<^C[G) 1^1^ — Z]cgC(G) ^(alnn)^ = (alnn)^. We obtain the following bound 
on err^: 

K a^ln^n l-K / ^dn^n, , , 

en.(/.K) < ^ _ p.^^ ^ p.^^ - = 0( — ) . (38) 

It is left to verify the conditions of Theorem 13. li Although we may assume that the statements 
of Theorem II . II hold with the respect to the approximate biased process, and \I{QK{t)) — y{t)\ = 
every < t < where y is the unique solution to ([T]), the 0{n ^/^) approximation 
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error will not be sufficient for proving that S{G) is approximated by z{G). The difficulty is in 
proving the Trend Hypothesis; dS{G)/dt is approximated by gi{I)+g2{I)S'^{G) for some continuous 
functions gi,g2, and using the y- approximation for I{G) implies an 0(n^^^^)S'^{G) error. Therefore, 
we apply Theorem 13.11 on both /(G) and S{G) (re-proving the result on I{G) for the interval 
[0,ti+i])- 

Keeping the notations of the theorem, define: 

Yi{G) = nI{G) = \Co{G)\ , 
Y2iG) = nS{G)= \C\\ 



fi{x,y,z) = 2( 



C6C(G) 

(.i-y)K 

y2 + (1 - y2)K 



hi^'V'^) = 2(^2^(1^^ 1) + 1) • 

The set T> will be a bounded connected open set containing the domain [0,tj+i] x [0,1] x [0,i], 

^ — r 

where i is the constant defined in (|3ip . Clearly, D contains (^, — ^^^^^) as long as T < tjj^in/2 
— T 

and S{Qk ) < ^, and in particular, for T = 0. 

Recall that the first coordinate of Y satisfies the Boundedness Hypothesis for a choice of /? = 
2 and every T < tj_|_in/2. Since every component of G is of size at most a Inn, (|32p implies 

■ — T+l — T 

that \Y2{Qk ) — Y2{Qk )| < 2(alnn)^, thus a choice of /3 = 2a^lnn verifies the Boundedness 
Hypothesis for every T < tj+in/2. 

The Trend Hypothesis was satisfied by Yi , fi in the proof of Theorem 11.11 using a choice of 
Ai = 1/n (/i is independent of z). By (f37|) and (p8]) . we have: 

IE (Y2(g^/^') - Y2{gK^)) - f2{^, ^^)\ = 2errM,K) = O(^) . 
V / n n n 

Therefore, a choice of Ai = 0(^^^) verifies the Trend Hypothesis. 

The Lipschitz condition is again satisfied by the fact that fi{x,y,z), i = 1,2, are clearly C°° 
for K >0. 

Notice that our definition of £, relying on Proposition 11.21 is such that the unique solution to 
the system of differential equations remains at least e-far from the boundary of the domain V for 
every < x < tj+i, hence a choice of cr = tj+i is valid. 

Altogether, a choice of A = implies that with probability 1 — o(exp(— in^/^)) there exists 

some constant C2 > such that the following holds for every < T < tj^in/2: 



n n n 



<C2n 



-1/4 
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where {y, z) is the (unique) solution to the equation system: 

^ = /i(^,y,^), ^ = /2(x,y,z) 
y(0) = z(0) = 1 

ScaUng the time to units of n/2 edges, we obtain that, with probabiUty 1 — o(exp(— ^n"*^/^)), the 
following holds for every < t < tj+i: 

\S{QK{t))-z*{t)\ = 0{n-^'^) , 

where z*{x) = z{x/2). Hence, z* is the unique solution to the equation: 

dz* _ldz _ K 2 (1 - K)y'^ 

Itx ~ 2dx ~ y2 _^ {l-y2)K^ ^ y2 _^ {I - y'i)K ' 

with the initial condition z*(0) = z(0) = 1, completing the proof of the lemma. ■ 
4.3 The susceptibility at the supercritical phase 

Proof of Theorem [1T3] part m Let Q = [max{i^:, ;^}] , and assume without loss of generality 
that e < tc- By Proposition 11.21 there exists some < < tc satisfying z{t~) = Lemma 14.31 

implies that, with high probability, 

S{GK{t-))>^{l-o{l))>^ , 

and in addition, Qxit') has a k,c component tail. Fix H = QK{t )- Setting i = and t = ^, 
Theorem 14.21 implies that Qi(-^)\i{ almost surely has a giant component. Hence, by Theorem 12.11 
(inequality dU), Gk{^)\h = QK{t~ + almost surely has a giant component, and in particular, 
QK{tc + e) has a giant component with high probability. 



4.4 Computer simulations and numerical results 

We have implemented an efficient simulator of the model Qk, which runs in a time complexity of 
0(n log n), and uses a data structure of size 0{n). Maintaining the sets of isolated and non- isolated 
vertices allows an 0(1) cost for randomizing the next edge at each step whenever one of its endpoints 
is isolated; the cost of selecting an edge between two non-isolated vertices is O(logn), where the set 
of existing edges is kept in a balanced tree, which is an efficient data-structure. Finally, maintaining 
the connected components in linked lists according to the Weighted-Union Heuristic (see, e.g., [8] 
p. 445) provides an average cost of O(logn) for uniting components. 

We conducted a series of simulations with n = 10^, and values of K ranging from to 10 in 
steps of 0.1. In order to assess the critical time tc{K), we considered the minimal time at which 
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Figure 2: Comparison of the numerical results for tc{K), and estimations of tc{K) according to 
computer simulations of the model. 

the largest component of Qk was of size an, where a = 0.01. These estimations of tc{K) were 
supported by tests of higher values of n, such as 2.5 • 10^ and 5 • 10^ vertices. We note that the 
constant a decreases with K (independently of n). Therefore, automatically assessing tc{K) for 
large values of K requires a different method, such as a threshold for the derivative of the size of 
the largest component, or for the derivative of the susceptibility. 

The simulated value of tc{K) in each setting was averaged over 10 tests conducted as mentioned 
above. Figure [2] shows an excellent agreement between these values and the numerical approxima- 
tions of tc{K), obtained by numerical solutions of the ODEs ([T]) and ([2|) by Mathematica. 

5 Analytical properties of the critical time 

In this section we prove Theorem II .41 regarding the analytical properties of the critical time tc{K). 
5.1 Asymptotical Behavior of tc{K) for K ^ 1 

Proof of Theorem II. 4[ The proof begins with the analysis of tc{K) as K ^ oo, whose leading 
order asymptotics is constructed by the method of matched asymptotics. 
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The differential equation ([T]) for y{x) 

r + (1 - r)^ 

with its initial condition y{0) = 1, suggests the transformation u = l — y,sou<^l near the origin, 
and satisfies 

— 7/ = — 

{I - + Kui2 - u) 

with the initial condition u{0) = 0. Hereafter ~ means equality of leading order terms. For n ^ 1 
we have 

, 1 + Ku 

whose solution is 

2Ku - ln(l + Ku) iTx. (41) 
For K^^ <^u the logarithmic term is smaller compared to the linear term, hence 

u(a;)«|, ioi K-^ <^ X <^ K-^ , (42) 

where e > is needed to ensure u <C 1. 

We rewrite the differential equation ([2]) for z{x) as 

z' = c(x)(z2-l) + l, (43) 

where 

+ (1 - y^)A 

We recognize (|43l) as a Riccati equation (see, e.g., |11J). However, it has no closed- form solution 
for a general function c(x), and we resort to asymptotic methods. 

We obtain an asymptotic approximation of z{x) by the method of matched asymptotics. As 
shown below, we first construct an asymptotic solution that is valid for <C x ^ ET"^/^, which 
is then matched to a different asymptotic solution that is valid for x 3> K~^/^ up to the blowup 
point. 

First, we use a linear approximation of ()43p near x = 0, by substituting v = z — 1 that satisfies 

v' = c(?j2 + 2v) + 1, (45) 

with the initial condition t)(0) = 0, so u is small near the origin. Therefore, as long as the nonlinear 
term satisfies cv"^ <C 1, it can be viewed as a perturbation to the linear equation 

v' -2cv = l + 0{cv'^). (46) 

Multiplying (06]) by its integrating factor we obtain 
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Equations (00]) and (jH]) imply 



Therefore, 



c(x) dx 

Hence, by equation (07|) we have 



Ku' 
^ l + Ku 

Ku' 



l + Ku 



dx = + Ku). 



v^{l + Kuf j — 



dx 



+ KuY ■ 

We calculate the integral in equation ([50]) by implying equation ([iO]) 



(48) 
(49) 

(50) 



dx 



[l + KuY 



l + 2Ku 



1 

'k 



l + Kuf l + Ku 

2 1 

+ 



dx 



[l + Kuy {l + KuY 



du 



+ 



Thus, 



l + Ku 2(1 + Kuy 2 

V PS ^u{2 + mu), 
which combined with equation (I42p for x S> i^^^ gives 



3 

v{x) PS -Kx^. 
8 

Equations (|44p and (|42p imply that c{x) is to leading order 



c{x) 



K 



-, for K ^ < a; < 



(51) 



(52) 



1 + 2Ku X ' 

Recall that the linear approximation for v{x) is valid whenever cv^ <C 1. Equations (|5ip and 
()52p give cv^ ~ -fC^x^, hence the linear approximation (|46p holds for x ^ il'"^/'^. Thus, the 
approximation (|5ip holds for i^"^ ^ x ^ 

Second, we choose an intermediate point xi such that K^^ <^ xi <^ A'^^/^, e.g., xi = 
for which c(xi) = ©(AT^/^), 7;(xi) = 0{K-^/'^), and z(xi)2 - 1 = 0{K-^/'^). Therefore, 
c(xi) (2;(xi)^ — l) = O (AT^/^) and we may approximate equation (jl3]) by neglecting the constant 
term 1 < K^l'^ 

z' ss c(x)(z2 - 1), (53) 



with the initial condition 



The blowup point tc satisfies 



The left hand side of (1551) is 



zi = z{xi) = 1 + v{xi) ?s 1 H — Kx\. 



dz 



Z^ - I 



c(x) dx. 



(54) 
(55) 



°° _dz 



_1, ( zi + l 

z^-l~ 2 \z^-lj' 
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whereas the right hand side of (j55|) is calculated by equation (j49j) and the approximation ([42 

1 + Kuc\ , / i + i/s:tc\ 



£c(x)<ix^ln(i±|j^)«ln 



1 + ii^Xl 



Therefore, 

zi + 1 ( 1 + iKt^ \ ^ tl 



zi-1 [l + ^Kxi 
Substituting and rearranging we obtain 



(56) 



2 ^2 + p<x{ 16 

^^'^'^^ iKxl ""SK- ^^^^ 
Note that the result is independent of the choice of the intermediate point xi, as expected. Hence 

tc ~ -^=, for K > 1, (58) 

is the leading order asymptotics of the blow up point. Throughout the asymptotic analysis we 
used merely the leading order terms in their validity regimes. Therefore, the higher order terms 
are guaranteed to be asymptotically smaller than the leading order term, that is, 

tc(i^) = ^(l+o(l)), forK»l. (59) 

We constructed numerical solutions of the ODEs ([I]) , (l2|) for various values of K using Mathemat- 
ica. Figure [3] shows the comparison of the numerical approximation for tc{K) and the asymptotic 
result A/^/m. 



5.2 Strict monotone decreasing behavior of tc{K) 

The continuity of the critical point tc{K) follows from the general continuous dependence of ODEs 
on their parameters, which is in our case a single parameter K. Thus, to complete the proof we 
must show that tc{K) is strictly monotone decreasing in K. 

We first show that c{x) is strict monotone increasing in K, i.e., for any < Ki < K2, the 
corresponding functions ci,C2 (where y is replaced by the corresponding solutions to ([T|)) satisfy 
ci(x) < C2{x) for every x > 0. 

Recall that by the properties of y, as stated in Proposition 13.31 c(0) = K and for every K ^ 0,1, 
c{x) either strictly increases or strictly decreases to 1 as x ^ cxd (depending on whether K > 1 
or K < 1). Hence, if Ki < 1 < K2 or Ki < 1 < K2, clearly ci(x) < C2{x) for every x > 0. 
Furthermore, if 1 < Ki < K2, the strict monotonicity of y in i^T implies that for every x > 0: 

_ 1 1 _ 
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Figure 3: Comparison of the numerical results for tc{K) obtained by numerical solutions of the 
ODEs ([H)-©, and the asymptotic approximation of Theorem II. 4[ Logarithmic scale was used in 
both axes. 

Assume therefore that < Ki < K2 < 1. Chain rule differentiation of (j44p . together with ([1]), 
gives: 

c'(x) = -2c' ((1 - y)c + 1) ■ (60) 



Rearranging equation (j44l) . we get: 

Combining equations ([60|) and (j6T]) . we obtain an equality of the form d (x) = g{c, K), where g{c, K) 
satisfies: 

dg , c^(l + c) l {l-c)K 

for all < ii' < 1 and x > 0. By the argument stated in the proof of Proposition 13. 3^ this implies 
that c is indeed strictly monotone increasing in K. 

Returning to equation (j43p . since z' = c{z^ — 1) + 1 and c is strictly monotone increasing in 
we conclude that z'[t) is monotone increasing in K for any t > (recall that, as stated in the proof 
of Proposition 11.21 •^'(i) > for any t > and hence z[t) > 1 for any t > 0). Hence, the above 
consideration implies that z is strictly monotone increasing in K. Therefore, tc{K) is monotone 
decreasing in K, that is, 

tciKi) > t,{K2), for Ki<K2. (63) 
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Yet, the fact that z is strictly monotone increasing in K does not guarantee by itself that tc{K) is 
strictly monotone decreasing in K. 

To this end, we consider solutions 2^1, 2:2 and ci,C2 corresponding to Ki < K2 and prove that 
tc{Ki) > tc{K2)- The functions ci and C2 are continuous on [0,oo) and satisfy ci < C2, therefore 
3e > such that 

C2(x) - ci(x) > e, Vx G [QMKi)\. (64) 
The function (whose choice is made clear below) 



-2. 1 
= 

m 

z 

satisfies lim2_»oo ^(-2) = 0. Since zi{x) ^ 00 as x ^ tc{Ki), there exists a point < t* < tc(-f^i)) 
such that 

zi(r) >max{l,i^fi} and .4(zi(r)) < ^i— -. (65) 

maxjl. All 

Let (rc) be the solution of the differential equation 

z = ci{x)z^ + ci(x)z, z(r) = zi{t*). (66) 

Recall that 21(2;) is the solution of 

z' = ci(x)z2 + l-ci(x), z{t*) = zi{t*). (67) 

Comparing and ([S7|) we note that ci(x) 2:^(2;) > 1 — ci(x) for all x > t*, because z*{x) is 
monotone increasing in x, ci(x) > min{l, Ai}, and the choice of t* in equation (j65p . indicated by 

ci(x)z^(x) > ci(x)2]'(t*) > min{l, Ai}max{l,A:f^} = 1 > l-ci(x). 

Hence, ^^(x) > 2:1 (x) for all x > t* and 

t,{Ki) > tliKi), (68) 

where t*(Ki) is the blowup point of zl{x). Similar considerations show that the blowup point 
t*{K2) of the solution, to the differential equation 

z' = C2(x)(z2_l), z{t*) = Z2{t*), (69) 

(remark: if t* > tc{K2) then tc{Ki) > tc{K2) and the proof is completed, so we assume t* < tc{K2)) 
satisfies 

tl{K2) > t,{K2), (70) 

because the positive term 1 was omitted in compared to (^5]) . The advantage of the differential 
equations ([66]) and ([69]) over the original equation ([l3j) is that their blowup points satisfy closed 
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form relations. Namely, 

= h"\.i.)d.. (72) 
Assume to the contrary that tc{Ki) = tc{K2). Then, by equations ()68p and (I70p it follows that 

t:(K2)>t:(^i). (73) 

Combining equations ([7T]) . ([72]) .([73 |) . ([M|) and the positivity of C2(a;), we obtain 



1 /z*(r) + i 



2^''V4M^J " C2{x)dx>e{tl{Ki)-t*)+ I ci{x)dx 



Employing equation (j7ip with ci(x) < max{l,i^i}, we find 



Equations ([75]) and ([77]) give 

4(t*)2 



In 
In 



zlit*) 

which contradicts ([65]) . Therefore, ic(-^i) > ic(-^2)- 



(74) 



Since In |^ is a monotonic decreasing function of z, and zl{t*) < Z2{t*), it follows that 
or equivalently 

zl{t*^^ 



^2;J(r) + l\ max|l,Ai} 



6 Concluding remarks and open problems 

As mentioned in |l3j, there is a well known relation between the "double jump" phenomenon in 
Qi and the critical threshold in percolation. Indeed, the susceptibility corresponds to the size of 
the component containing the origin, and the appearing of giant component is analogous to the 
appearing of an infinite cluster. The subcritical and supercritical stages in the Qi correspond to 
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those in percolation. However, it is unclear how to define a percolation process which is analogous 
to the biased process. 

We conclude by describing several natural ways to extend the family of biased processes and 
their analysis. 

Consider processes where the probability of adding an edge at each step is biased according to 
the sizes of the components of its endpoints. For instance, it is possible to assign the weight Kij to 
a missing edge {u,v), such that \Cu\ = i,\Cu\ = j for i,j < C, where C > is some large constant, 
and some weight Kq otherwise. It seems plausible that the above methods can be applied in this 
scenario as-well. 

The biased process Qk was defined as a graph process on the complete graph on n vertices, 
Kn, in other words, at each step we considered the entire set of missing edges, and the process 
ended when obtaining Kn- It is possible to define an analogous process on a different underlying 
graph H = (V,E), where the edge pool at each step contains only edges belonging to E, and the 
process ends once H is obtained. For instance, it is possible to consider the biased process on Qn, 
the n-th dimensional cube, or L2{n), the n xn grid. In these cases, it may be interesting to study 
the critical time for satisfying other monotone increasing properties, such as left-right crossing in 
L2{n). It seems plausible that this crossing probability is monotone decreasing in K. 

Acknowledgement We would like to thank Itai Bcnjamini, who proposed the model of the biased 
process and provided useful insights to the problem. 
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